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ABSTRACT 

Context. Age, metallicity and spatial distribution of globular clusters (GCs) provide a powerful tool to reconstruct major star- 
formation episodes in galaxies. IKN is a faint dwarf spheroidal (dSph) in the M 81 group of galaxies. It contains five old GCs, 
which makes it the galaxy with the highest known specific frequency (S^v = 126). 

Aims. We estimate the photometric age, metallicity and spatial distribution of the poorly studied IKN GCs. We search SDSS for GC 
candidates beyond the HST/ACS field of view, which covers half of IKN. 

Methods. To break the age-metallicity degeneracy in the V-I colour we use WHT/LIRIS Ks -band photometry and derive photometric 
ages and metallicities by comparison with SSP models in the V, /, Ks colour space. 

Results. IKN GCs’ VIKs colours are consistent with old ages (> 8 Gyr) and a metallicity distribution with a higher mean than typical 
for such a dSph ([Fe/H] =: -1.4 /q 2 dex). Their photometric masses range (0.5 < Aloe < 4 x 10^ Mq) implies an unusually high mass 
ratio between GCs and field stars, of 10.6%. Mixture model analysis of the RGB field stars’ metallicity suggests that 72% of the 
stars may have formed together with the GCs. Using the most massive GC-SFR relation we calculate a SFR of ~ lOA/g/yr during 
its formation epoch. We note that the more massive GCs are closer to the galaxy photometric centre. IKN GCs also appear spatially 
aligned along a line close to the IKN major-axis and nearly orthogonal to the plane of spatial distribution of galaxies in the M 81 
group. We identify one new IKN GC candidate based on colour and PSF analysis of the SDSS data. 

Conclusions. The evidence towards i) broad and high metallicity distribution of the field IKN RGB stars and its GCs, ii) high fraction 
and Hi) spatial alignment of IKN GCs, supports a scenario for tidally triggered complex IKN’s SFH in the context of interactions with 
galaxies in the M 81 group. 

Key words, globular clusters: general - galaxies: individual (IKN) 


1. Introduction 

Globular clusters (GCs) are compact stellar systems which typ¬ 
ically form on a short time-scale, and as such, their stellar pop¬ 
ulations are composed of stars with similar age and metallicity. 
Therefore, the GCs’ integrated light can be closely approximated 
by a single stellar population (SSP) with a given age and metal¬ 
licity. This makes GCs unique tracers of past major star forma¬ 
tion episodes in a galaxy’s star formation history (SFH). They 
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IKN is a faint, low surface brightness dwarf spheroidal 
(dSph) galaxy {My - -11.5 mag, i.e. Ly — 3.4 x 10 ^Lq). It 
is located in the outskirts of the nearby M81 group of galaxies 
(m-/W = 27.79+ 0.02 mag, 3.61 Mpc, |Dalcanton et al.|2009| ), and 
at about 82 kpc projected distance from M 81. Recent analysis of 
Hubble Space Telescope (HST) imaging data with the Advanced 
Camera for Surveys (ACS) revealed that it has five GCs whose 
V-I colours and half-light radii are consistent with being old 
and metal-poor ( Georgiev et al.||2()09b| l. The brightest of these 
GCs, IKN-05, was indeed spectroscopically confirmed to be a 
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Such high S n therefore poses the puzzling question of what 
mechanisms could be responsible for the formation of the unusu¬ 
ally large number of GCs that make up 13% of IKN’s luminos¬ 
ity? A plausible explanation could be offered by a complex for¬ 
mation history of the IKN, triggered by tidal interactions (tidal 
formation?) with more massive galaxies in the M 81 group. Such 
a scenario is supported by the finding that IKN field stars have a 
broad metallicity distribution, which peaks at a too high metal¬ 
licity for its luminosity, [Fe/H] = -1.08 ±0.16 dex (Fianou et al. 
2010| l. Dwarf galaxies with similarly high metallicity for their 
luminosity have been also observed in the M 81 group and were 
also suggested to have formed during galaxy interactions in the 
group from tidal debris (e.g. Croxall et al.||2009 l. Here we aim 
to test this scenario by studying the integrated properties of the 
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IKN GCs. In addition, the high fraction of stars in GCs carries 
important information and implications for questions like how 
many stars formed in clusters in a given star formation epoch 
and how many clusters remain bound, i.e. implications for star 
cluster dissolution (e.g. Lada & Lada|2003] Bastian|2008[ God- 


dard et al.|2010[[Kruijssen|2012 Larsen et al.|2014a|r 


In this study we aim to improve the age and metallicity 
information for all GCs in the IKN dwarf spheroidal (dSph), 
by combining optical and near-infrared (NIR) photometry. We 
use existing optical (HST/ACS, |Georgiev et al.|2009bl l and new 
near-infrared K^-band data from LIRIS on the William Herschel 
Telescope (WHT). The latter is very important to break the age 
metallicity degeneracy in the optical. This problem can be par¬ 
tially resolved with a combination of optical and NIR colour in- 
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therein). Additional complications in the optical domain come 
from the hot evolved stars (HB, blue stragglers) which contribute 
significantly to the integrated light, whereas in the NIR and for 
ages > 2 Gyr, the red giant stars completely dominate the total 
energy emission. In effect, NIR colours such as I-K or J-K mea¬ 
sure the effective temperature of red giant branch stars, which 
in turn is very sensitive to metallicity, without a strong age de¬ 
pendence ( |Worthey et al.|[T994| l. An advantage of the optical- 
NIR technique, when compared to spectroscopy, is that for lesser 
amount of observing time, the properties of a statistically repre¬ 
sentative sample of GCs can be estimated. In galaxies beyond the 
Local Group resolving stars is very difficult with current ground- 
based telescopes. However, GCs are detectable and can be stud¬ 
ied as far as the Centaurus, Hydra I and Coma galaxy clusters 
(ss 100 Mpc away) (e.g. |Peng et al.|2011| l, thus providing a re¬ 
liable and powerful tool to study galaxy formation in the local 
universe. In addition, knowing about the number of GCs helps us 
shed light on galaxy formation theories. Therefore, we perform 
a search for new GC candidates in SDSS around IKN. The GCs 
are expected to differ from stellar point spread functions (PSFs) 
due to their more extended nature (effective radius of about 3 pc) 
and thus would appear marginally resolved even in ground based 
imaging, such as the ones from SDSS. 

The outline of the paper is as follows: in Section]^ we de¬ 
scribe the NIR-Optical data set and the reduction procedures. 
In the sections in §[^we analyze and discuss the implications 
for the formation of IKN based on the photometric properties 
of the IKN GCs from derived photometric age, metallicity and 
mass (§|3.1|l, SFR (§ 3.2 1 , mixture models for the metalicity dis¬ 


tribution of IKN RGB stars (§ |3.3| l and GCs’ spatial distributions 
(§|3.4|and §|3.5[). In §[^we summarize our results. 


2. Observations and data reduction 

2.1. Observations 

The IKN dSph was observed in the K^-band (2.15/im) with 
LIRIS (Long-slit Intermediate Resolution Infrared Spectro¬ 
graph), mounted on the 4.2 m William Herschel Telescope 
(WHT), during two consecutive nights (15 and 16 March 2009). 
LIRIS provides a field of view of 4.27 x 4.27 arcmin^, thus cov¬ 
ering the entire galaxy and its immediate surroundings. A plate 
scale of 0.25 "/pixel (^ 4.5pc/pixel at distance of 3.61 Mpc to 
IKN) means that its globular clusters will be observed as unre¬ 
solved extended sources. 

The Kj-band observations (see [Chies-Santos et al.|[2011b| l 
were performed in a dither pattern of five consecutive exposures 

Article number, page 2 of 1 11 1 


at five different positions in a cross-like arrangement with offsets 
between 10" and 20". The exposure time of a single image (DIT) 
was kept short (15 seconds) to deal more easily with the high sky 
background level in the near-IR. This observing strategy facili¬ 
tated the sky subtraction reduction step (see Section 2.2. l| i. 

The final exposure time of all combined images sums up to 
5505 seconds, which allowed us to perform photometry on the 
faintest IKN GC (V - 21.2 mag, i.e. My — -6.65+0.06 mag) for 
a distance modulus adopted here of {m-M) - 27.79 + 0.03 mag. 
The observation log is summarized in Table[T] A total of 6 photo¬ 
metric standard stars for photometric calibration were observed 
in each night. 


Table 1 Observations log and night conditions. 


Night 

FWHM 

N 

Total exposure time 

yyyy-mm-dd 

" 

images 

seconds 

2009-03-15 

1.0 

76 

1140 

2009-03-16 

=s0.9 

291 

4365 


The V and I photometry of the IKN GCs ( |Georgiev et al.| 
2009b|l used in this study is based on the HST/ACS archival data 


program SNAP-9771 (PI: I. Karachentsev). The non-dithered 
2x 600 s 7^6061^ and 2 x 450 s F814W exposures were designed 
to reach the Tip of the Red Giant Branch (TRGB) for distance 
measurement ( [Karachentsev et al.|2007| l. 

To extend the search for GCs around IKN in regions not 
covered by the HST/ACS observations, we used the Slo an Dig 
ital Sky Survey (SDSS) ugriz DRIO archival image^ (Abaza 


[jian et ak 2009 | l. SDSS is a large scale survey with the 2.5m 
Apache Point Observatory telescope covering a quarter of the 
sky. SDSS reaches a resolution of 0.396"/pixel with an r-band 
seeing of 1.55" and limiting magnitudes sufficiently deep to en¬ 
able a search for GC candidates. The pixel scale corresponds to a 
projected spatial resolution at the IKN distance of 7 pc/pixel and 
the corresponding seeing value of 1.55" to a Gaussian FWHM 
of 13.5 pc. The ugriz data allows us to extend the wavelength 
coverage for the already known clusters and search for new can¬ 
didates (see details in Section |2. 2. 3| l. All magnitudes are in the 
Vega system. 


2.2. Data reduction and photometry 
2.2.1. LIRIS data reduction 

The NIR imaging reduction was performed within the LIRIS 
IRAF package LIRISDR]^ For each night the data was first 
corrected for the pixel mapping anomaly with the Icpixmap 
IRAF task, an effect produced by the reconstruction of the two- 
dimensional array which misplaces some of the pixels, espe¬ 
cially in the lower left quadrant. 

The images were flat-held corrected and sky subtracted us¬ 
ing the Idedither LIRISDR task for each night. One of the first 
corrections performed is the flat-fielding, which must correct not 
only for large scale gradients in the image but also for pixel-to- 
pixel gain variations. The dome flat fields were obtained by sub¬ 
tracting dark frames from bright frames of the same exposure 
time in order to suppress the possible thermal contamination of 
the telescope environment and then coadded and normalized for 

* http://www.sdss.org/drlO/ 

^ Package written by Jose Acosta Pulido to man¬ 
age data produced by the LIRIS instrument 
www.iac.es/galeria/jap/lirisdr/LIRIS_DATA_REDUCTION.html 
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Fig. 1 SDSS gri colour composite mosaic of the IKN dSph field and its five known GCs. These are indicated with squares in 
the colour image and shown magnified on the right from the HST/ACS image data. The labels indicate the field of view for the 
archival HST/ACS FhOhlT, FSITIT and our WHT/LIRIS K^-band image. The light-colour stellar density contour map is adapted 
from Figure 6 in Lianou et al. ( |2010| l. The ellipse indicates the HyperLEDA values for the IKN diameter (2.69' =2.91kpc) and 
ellipticity (0.18). For reference, with the arrows we indicate the direction to the group dominant galaxies M 81 and M 82. 


each night. The sky subtraction step is the most important step in 
the data reduction of NIR data as the brightness and structure of 
the NIR sky varies on a timescale of minutes and its contribution 
must be individually subtracted for each pixel. The general idea 
is to obtain the median frame for a set of images, discarding the 
values of the pixels that have received flux from celestial objects 
and then subtract the median frame from the science frames. The 
Idedither procedure creates sky frames from sets of five consec¬ 
utive (dithered) exposures and generates a sky subtracted image 
for each frame before coaddition. 

We aligned all non-registered, but sky-subtracted output im¬ 
ages from Idedither using a custom IRAF procedure based on 
geomap, geotran and imcombine IRAF tasks. We used a sigma¬ 
clipping coaddition algorithm, with a statistics region used to 
estimate the images zero offsets and the images are scaled by 
their exposure time. Figure[2 shows the final overlap region of 
the LIRIS coadded image with the HST and SDSS DRIO data. 
The standard stars data were reduced using the same IRAF and 
LIRISDR procedures as for the IKN images. The more photo¬ 
metric night was March 16, to which all images were subse¬ 
quently aligned and registered. The total exposure time in the 
final coadded image is 5505s. 


2.2.2. Ks-band photometry of IKN globular clusters 

Aperture magnitudes of the clusters, standard and held stars 
were measured with tasks in the IRAF DAOPHOT package. We 
performed a curve of growth analysis with r -2, 3,4, ...,40 pix¬ 
els aperture radius for the GCs and high signal to noise isolated 
held stars in the images, subsequently choosing a photometry ra¬ 
dius of 12 pixels. This is about three times the average full width 
half maximum of stars and unresolved sources in the final co¬ 
added image. The curve of growth analysis from isolated stars 
allowed us to derive an aperture correction of 0.04 magnitudes. 
PSF photometry was not possible due to the low number of good 
PSF stars necessary to build a reliable PSF model. 

The photometry is calibrated assuming a transformation 
equation of the form below, where the coefficients are derived 
from a least squares fit to the magnitudes of the photometric 
standards; 


Ks ^ks + Zk+ Pk^ Xk^ff + c X (V - ks). (1) 

Ks and ks are the standard and instrumental magnitudes, pk = 



Article number, page 3 of 1 11 1 
















A&A proofs: manuscript no. main 


( |1990| l. Due to the low number of /T-band standards with reli¬ 
able published V-band magnitudes, the value of the colour term 


c was adopted to be 0.02, as found in earlier studies (Goudfrooij 


|et al.|200T||Georgiev et al.|2012| l. The zero point was evaluated 
from the fit to the standard stars to be 23.01 + 0.02mag, not too 
different from the value found by ( Chies-Santos et al.|2011b I of 
23.07 ± 0.01 mag for another set of observations. 

The final calibrated ^f^-band magnitudes of the clusters 
were estimated by using their V-band magnitude published in 
Georgiev et al. ( 2009b| l, iterating Equation[T]until the Ks magni¬ 
tude converged. The V and magnitudes were also corrected 
for foreground galactic extinction with values for the total ab¬ 
sorption of Ay = 0.22 mag, A^s = 0.02 mag taken based on 
[Schlafly & Finkbeiner| ( |20lT I re-calibration of the |Schlegel et al.| 
( [l998| l dust maps assuming the |Fitzpatrick ( T999] l reddening law 
with f?i/ = 3.1. 

The foreground reddening corrected IKN GCs photometry is 
tabulated in Table|4] 


2.2.3. Search for GC candidates in SDSS DR10 

Since HST/ACS field of view covers half of IKN, we utilize the 
wide spatial and spectral coverage of SDSS to search for GC 
candidates. IKN dSph is not entirely covered by the HST/ACS 
held of view and the stellar density map ( Lianou et al.|201()| ) in¬ 
dicates that the center of the IKN held stars is approximately at 
the location of the bright foreground star near the edge of the 
held of view (see Fig.[T] also § |3.4| i. Therefore, given the appar¬ 
ent distribution of GCs on one side of IKN additional GCs may 
be located outside the HST/ACS image. Searching for GC can¬ 
didates beyond the optical extent of IKN is also motivated by the 
discovery of very distant GCs around the dwarf irTegular galaxy 
NGC 6822 located out to 11 kpc from the galaxy centre ( [Hwang 
et al.|20lT] l. These GCs have similar luminosities, sizes and spa¬ 
tial distributions as the known IKN GCs (see §|^and § |3.5] l. 

To search for GC candidates in IKN we make use of their 
nature as extended objects. HST/ACS derived half-light radii for 
the known GCs are in the range 1.96 + 0.16 pc to 14.81 + 0.83 pc 
( [Georgiev et al.|2009b| l. The r-band seeing of the SDSS images 
used is 1.55", which corresponds to a projected spatial scale at 
the IKN distance of 27 pc. This implies that at least the larger 
half-light radii GCs will have a PSF which differs from a point- 
source PSF. 

Although the model and PSF magnitudes in the SDSS pho¬ 
tometric catalogues are quite reliable, we perform our own PSF 
photometry on archival SDSS images in order to evaluate the 
difference between stellar PSF and GC profiles and search for 
new cluster candidates. The additional reason to perform our 
own PSF analysis is because IKN falls on two SDSS stripes, 
which requires a separate PSF analysis. We hnd the difference 
between the PSFs of the two stripes to be well within the mea¬ 
surement uncertainty. We selected isolated, high S/N stars to 
build a PSF model and performed a PSF photometry on all de¬ 
tected objects with the IRAF’s DAOPHOT package procedures 
ipsf and alltstar). Instrumental PSF magnitudes were calibrated 
to standard SDSS magnitudes by using the SDSS PSF magni¬ 
tudes of the same stars. The final SDSS photometry of the IKN 
GCs is tabulated in Table[4l 

Due to the fact that the clusters were not detected or with 
too low S/N in the M-band images, this band was not used in 
the subsequent analysis. We used the sharpness value returned 
by allstar during the PSF fitting, to search for objects more 
extended than the stellar PSF with similar sharp-walue as the 
known GCs. The allstar sharpness parameter is defined as the 
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Fig. 2 Selectinon of GC candidates based on sharpness value. 
Shown is a g-band example of the sharpness value versus mag¬ 
nitude for all sources with PSF photometry (small crosses), the 
known IKN GCs (large solid, blue circles) and the SDSS GC 
candidates (open circles). The PSF stars are shown with small 
solid symbols (dark circles). 


difference between the squares of the width of the PSF model 
and the width of the object. That is, sources with light profile 
similar to the stellar PSF will have sharp-value close to 0, while 
more extended sources will have a sharp-value greater than 0. 
The extended nature of the candidates can be seen in Figure]^ 
where the sharpness values for PSF stars, GCs and all other ob¬ 
jects are plotted as a function of their g-band magnitude. The 
second selection criteria was based on colours, with the candi¬ 
dates having to be in the immediate vicinity of the GC’s in the 
griz colour space (within 2cr from the average GC colours). 

One of the new GC candidates is very close to the bright 
foreground star (BD69-I-557), which prevents reliable photome¬ 
try. Therefore we do not consider it for further analysis. 


3. Analysis and Discussion 

3.1. GCs’ colours, photometric age, mass and metatlicity 


To gain insight on the formation epochs of the IKN GCs and 
the evolutionary path of this faint dSph, we discuss the stellar 
population age and metallicity of the IKN GCs based on their 
optical-NIR integrated colour distributions. As mentioned ear¬ 
lier, we take advantage of the good age and metallicity resolution 
provided by the optical-NIR colour-colour indices. 

Figure[^shows the iV-I)o vs {I-Ks)o colour-colour distribu¬ 
tion of all GCs in the IKN dSph, compared to Bruzual & Chariot 
P003[ l SSP models for a canonical ( [Kroupa[2001[ l IMF, and a 
range in age and metallicity as indicated by the labels. The dis¬ 
tribution of the GCs’ colours in Figure|^suggests that all of them 
are likely old and metal-poor [M/H]< -0.7 dex. This is con¬ 
sistent with the metal-poor and old age inferi'ed for the bright¬ 
est IKN GC-5 from recent spectroscopic analysis ( [Larsen et al. 
[2014a[ ). However, given the photometric uncertainties, ages as 
young as 2-5 Gyr cannot be ruled out. In general, although with 
lesser certainty, the photometric metallicity spread of the IKN 
GCs is consistent with that of the IKN held stars, [Fe/H]< -1.0 
dex ( Lianou et al.|2010| . We have also compared the IKN GCs 
with the GC systems of the luminous elliptical galaxies M87 and 
M60 from[Chies-Santos et al.[([201 laji in a gzk colour-colour di- 
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Fig. 3 Foreground redenning corrected optical - NIR colour- 
colour plot, (y-/)o V {I-Ks)o of the IKN globular clusters (solid 
labeled circles). The cluster colours are compared to |Bruzual &| 
Chariot p003| l SSP model isometallicity tracks (M/H labeled 


lines) and isochrones (age labeled lines). The top left arrows 
show the direction of increasing age and metallicity and the ar¬ 
row to the right shows a reddening vector of Ay = 0.5 mag. 


agram. IKN GCs 4 and 5 fall on the same loci as M87/M60 GCs 
in the gzk parameter space, implying old ages. The other IKN 
GCs do not have gzk photometry. 


To derive photometric ages and metallicities for the GCs in 
our sample we perform a minimization interpolation between the 
age and metallicity in the |Bruzual & Charlot| ( [2003] ) SSP model 
grid. These are estimated with the average value between the two 
closest values weighted by the distance to the respective model 
and the photometric uncertainty. This method enables us to de¬ 
rive upper and lower errors for the photometric ages and metal¬ 
licities. The validity of this approach has been discussed in de¬ 
tail for the M 31 and MW GC colours with IBruzual & CharloTl 


( 2003|l SSP tracks for old ages (e.g. |Puzia et al.|2002[ Georgiev 


|et al.|2012| ). GCs falling outside of the model grid have age and 
metallicity derived from the nearest age and metallicity model 
tracks. A comparison of this technique for different SSP model 
tracks is presented in [Georgiev et al.|P012| l. Errors due to the 
choice of the SSP model are included in our measured values 
and amount to about 10% at given age and metallicity. Thanks 
to the combination of optical-NIR magnitudes, uncertainty in de¬ 
rived cluster mass due to stochastic effec ts are below the 1-5% 
level for old and massive GCs, > 10^Mq (Fouesneau & Lan 9 on 


|2010[|Fouesneau et al.|2012| l. In Table|^the measured values for 
photometric age, metallicity and mass of the IKN GCs are tabu 
lated. 


Table 2 Photometric metallicity, age and mass of IKN GCs. 


ID 

Z 

[Fe/H] 

Age 

Mass 


[Zo] 

[dex] 

[Gyr] 

[10^ Mo] 

IKN-01 

IKN-02 

IKN-03 

IKN-04 

IKN-05 

u.uz.h-_oqi5 

o.o5i!°:«« 

0.078!;!;“^ 

0 069^" 

1 o'5’:2+0.235 
^•o-’^_0.423 

-L506+^-3| 
-1 32 

i.jz,i_0 442 

-1 

lO / J_oggO 

14 77+1-14 
' -1.30 

15.51+“« 

13.21+™ 

14.19!«“ 

13.80+^“ 

0 771 +0 069 

-Q-Q25 

1 

'^-0.299 

3.824/{;|’ 


To convert from colour-magnitude to metallicity-mass distri¬ 
butions, we used the SSP model MjLv values based on the GC’s 
colours to convert from luminosities to masses (for the adopted 
distance to IKN, see §[^. The resulting mass-metallicity distri¬ 
bution is shown in Figure]^ The absolute values for the metallic¬ 
ities of the IKN GCs range from -1.6 dex to -1.1 dex and are rela¬ 
tively metal-poor, and are also consistent with metallicity values 
of the IKN foreground stars (Lianou et al.|201())l. 



[M/H] [dex] 


Fig. 4 Mass-metallicity (top) and age-metallicity (bottom) distri¬ 
butions of the IKN GCs (solid labeled circles) derived from their 
VIK magnitudes. 

Figurej^ shows that the photometric GC ages are older than 
about lOGyr and with a narrow spread. This suggests that the 
early formation of all clusters occurred within a relatively short 
time scale, which is typically the case for the old GCs in other 
galaxies, such as in the Milky Way or the old GCs in the EMC 
(see more detailed discussion of IKN RGB and GCs’ metallicity 
distributions in §|3.3|). 


3.2. IKN peak SFR from its most massive GC 

It has been demonstrated that there is an observational relation 
between the galaxy star formation rate (SFR) and the luminosity 
of its most massive cluster at a given epoch ([Larsen & Richtler 


2000 Weidner et al.|2004| Maschberger & Kroupa|2007[ Bastian 


2008 1 . Here we derive the IKN peak SFR from the photometric 
mass of the IKN most massive GC identically to [Georgiev et al.[ 
( 2012| l. Briefly, this included a correction for globular cluster 
mass loss due to stellar evolution using the [Bruzual & CharloT| 
( [2003 [ l SSP MjLy ratios at 10 Myr. We have not accounted for 
mass loss from tidal shocking, which however, could be expected 
to be small for such a low mass dwarf. Because the suggested 
age spread among the IKN GCs is very small (cf Fig.|^, the 
IKN peak SFR is calculated from its most massive cluster. For 
reference, we show in Figure]^ the age against the SFR derived 
for each GC. The upper limit estimate for the peak star forma¬ 
tion is of 10 Mo/year and occurred around 14 Gyr ago. The most 
massive GC provides a lower limit of the SFR at 14.29]^['2} Gyr, 
(the errors are propagated from the GCs mass uncertainties). The 
duration of the star formation burst cannot be estimated with 
high precision without good knowledge the GCs ages. However, 
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Fig. 5 Star formation rate calculated for each IKN GC against 
cluster age. The most massive GC (GC-5) yields a SFR of 
10 Mq! year at the time of GCs formation. 



given the low mass of the IKN and assuming a constant SFR, 
we can constrain the maximum length of the burst to a only a 
few (ten) million years. This indicates that IKN went through a 
bursty star formation episode at the epoch of the GCs formation. 

Due to the very small number of globular clusters, one must 
keep in mind that statistical effects in the GC mass distribution 
might play a significant role in the SFR that we find for IKN. 

The new GC candidates are either outside the HST and WHT 
field of view or within the region rendered unusable for analysis 
by the bright star. Thus, the candidates could not be confirmed 
visually on the HST image and that the NIR data is not available 
for breaking the age-metallicity degeneracy. 


3.3. Metallicity distributions of GCs and fieid RGB stars - 
impiications for IKN SFH 

The comparison between the metallicity distributions of the GCs 
and IKN’s field stellar population and using the knowledge of the 
galaxy SFR, can provide important insights about how much of 
the initial star formation took place in star clusters. 

In order to be able to directly compare the metallicity distri¬ 
bution function (MDFs) of the IKN GCs with that of the RGB 
field stars we convert their total metal content ([M/H], derived in 
§ 3.1 from SSP models) into [Fe/H]. For that we use the |Salan^ 
|et al'l ( |1993| l relation between [M/H] and [Fe/H] (their eq. 3), 
where the factor fa - 10“ with a — 0.3 dex, accounts for the 
typical high a-element abundance (i.e. low Fe content) observed 
in GCs (cf eq. 1 in [Ferraro et al.|1999| l. The converted [M/H] to 
[Fe/H] metallicity values for the IKN GCs are given in Table]^ 

Figure]^ shows the direct comparison between the MDFs of 
RGB field stars and IKN GCs. For viewing purposes only, the 
y-axis values of the GCs are obtained by multiply ing their pho 


Lianou et al. 


tometric mass, Mac, (cf Table|2| by 1.5 x 10“ 
kindly provided us with the photometry of the IKN field RGB 
stars, which MDF is shown with a solid line histogram in Fig¬ 
ure]^ [Lianou et al.j ( [2010[ ) have already reported that IKN has 
a broad MDF which peaks at higher [Fe/H] than for such low- 
luminosity dwarf galaxy. This suggests that IKN has likely ex¬ 
perienced a complex SFH with probably more than one burst of 
star formation. 

In order to assess the significance and the likelihood for a 
multimodal as opposed to a unimodal MDF and compare the 
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Fig. 6 Top: Metallicity distribution of IKN field RGB stars 
(solid, red line histogram) and IKN GCs (labeled solid circles). 
For viewing purposes, the vertical position of the GCs is scaled 
according to their mass. Thick solid (black) and dashed (grey) 
curves are the scaled probability density estimate (p*_rgb) and 
its Icr estimated from the model PDFs (thin gray lines, details 
in § [3.3| l. Bottom: Mixture model of the IKN field RGB stars 
metallicity distribution. The best fit maximum likelihood mix¬ 
ture model probability density distribution is shown with the 
solid blue curve and its components by filled curve Gaussians. 
For comparison we show again the metallicity distribution of 
IKN field RGB stars (solid red line histogram). The percent frac¬ 
tion of each fitted Gaussian is given in the legend and Table[^ 


fraction of stars in GCs, we performed a more detailed statisti¬ 
cal analysis with packages in r[^ First, we estimated the [Fe/H] 
probability density function (PDF) of the RGB stars using an 
Epanechnikov kernel density estimator. We used this PDF as a 
prior and the [Fe/H] uncertainty of the RGBs to generate a large 
number of models (1000). The different model realizations are 
shown with thin (grey) lines in Figure]^ The joint PDF of all 
the models, as a function of [Fe/H], gave the final PDF of the 
RGB stars [Fe/H] distribution. We show this in Figure]^ with 
solid (blue) curve and its corresponding one sigma intervals are 
represented by dashed (grey) lines. This test confirms that there 
is a dip in the IKN field RGB’s MDF at the Icr level. 

In order to further test the significance of the Icr dip in the 
IKN MDF, we perform an additional mixture model statisti¬ 
cal analysis to test whether the MDF is consistent with one or 
multiple underlying distributions. This particular test will allow 
us to assess whether the observed MDF can be represented by 
the combination of more than one Gaussian distribution. This 
will allow to derive the fraction of the stellar population which 
falls into the same metallicity range as the old globular clus¬ 
ters, i.e. within the same SF epoch. For that we used routines 
in the mixdist and mixtools R packages with which we fit a 


^ R is a free software environment for statistical computing. 
http://www.r-project.org/ 
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Table 3 Coefficients of the different Gaussian mixture model 
components. Column (1) gives the fraction of RGB stars in each 
of the components; column (2) and (3) give the Gaussians ju and 
Q- with their uncertainties. 


Population fraction 
[%] 

(1) 

P 

[dex] 

(2) 

cr 

[dex] 

(3) 

9.45 + 5.1 

-2.11+0.084 

0.1497 + 0.039 

36.09 + 16.3 

-1.57 + 0.129 

0.2199 + 0.067 

26.55 + 23.4 

-1.27 + 0.041 

0.1276 + 0.077 

27.91+9.0 

-1.00 + 0.023 

0.0837 + 0.011 


mixture Gaussian distribution models to the data by the method 
of maximum likelihood minimization. We performed a mixture 
models with one, two, three, four, five and six Gaussian com¬ 
ponents. The four Gaussian model was the preferred solution 
with considerably the best statistics (smallest highest 

p-value) of all models. Models with higher number of Gaus¬ 
sians give a marginally better x^^ but at the same time give the 
worst p-value. This is because all components are simultane¬ 
ously fit and likelihoods maximized. We can thus conclude from 
our statistical test that more than one component is needed in or¬ 
der to reproduce the MDF of the IKN RGB stars. We also note 
that only the four component model gave a good representation 
of the metal poor peak at [Fe/H] ^ -2.2 dex. It is beyond the 
scope of this paper to further investigate exactly how many com¬ 
ponents, i.e. episodes of star formation could have been respon¬ 
sible for the present day IKN field RGB’s MDF. The values for 
each of the four fitted components are given in Tablej^ 

Our mixture model fitting suggests that nearly a third (about 
28%) of the IKN field stars are in the most metal-rich compo¬ 
nent, while the rest 72% of the metal poorer RGB stars could 
have formed in more than one SF event. The uncertainties in 
the parameters of the three metal-poor components (see values 
in Table[^ prevent us from a more firm and detailed discussion 
whether and how many episodes of SF IKN could have went 
through. In addition, whether a Gaussian or another function 
would be a better description is also beyond the scope of this 
paper. For the following discussion we will consider that up to 
72% of the IKN star formation occurred simultaneously with the 
epoch at which the GCs formed. This stems from the distribu¬ 
tion of the IKN GCs in Figurej^a) compared to the field stellar 
component of IKN, which fall in the metal poorer component. 


We can first calculate the fraction of mass in clusters to to¬ 
tal galaxy stellar mass. Summing up cluster masses in Table]^ 
gives Adcc.tot = 8.14 x 10 ^Mq. For a total IKN luminosity of 
My - -11.2mag and a MjLy - 3 for [Fe/H] = -l.Odex and 
an age of 12Gyr, the IKN mass is thus A1* ik:n = 7.7 x 10 ^Mq. 
Therefore we obtain a mass fraction of /m.gc = 0.106, i.e. 10.6% 
of the IKN stellar mass is in its old GCs. This is unusually high 
compared to the fraction of light or mass in GCs in more mas¬ 
sive early- and late-type galaxies (e.g.|Peng et al.|2008] Georgiev 


|et al.|[20T0{ [Harris et al.|2013| l. If we consider our results from 
the mixture models above, i.e. that probably about 72% of the 
stars in IKN have formed in the same star formation episode as 
the old GCs, we will arrive at an even larger fraction of about 
14%. 


It is interesting to ask how the fraction fM,cc relates to the 
actual fraction of stars that formed in clusters, F = CFR/SFR, 
introduced by Bastian ( 2008| l. We can use the SFR =10 Mo/yr 
value obtained in § |3.2| and calculate the surface SFR, Ssfr- We 
adopt a radius equal to the optical extent of IKN, r = 1.5 kpc. 







0.5 1.0 1.5 2.0 2.5 3.0 

X (Kpc) 


Fig. 7 Two dimensional number density and histogram projec¬ 
tion distributions of the IKN field stars. The histograms have a 
bin width of 0.345 kpc and are fitted with a Gaussian function 
shown with solid lines. The outer countour line in the 2D plot 
indicates a stellar number density of 17.5 ★ /arcsec^ (16 * /pc^). 
The white star indicates our best fit position of the IKN centre 
with an uncertainty of + 150 pc. 


Thus, SsFR = IAMqI yr/kpc. This value for Ssfr suggests that 
F should be in the range F = 10 - 30% based on the F - Ssfr 
relation recently investigated for a wide range of environments 


(e.g. Bastian 2008| Goddard et al.||20T0{ Kruijssen|2012[ Mora 

jet al.||2014| ). This further confirms that IKN had a very intense 

SFH that lead to high CFRs and to the very high present day 
fraction of old GCs. 


3.4. Spatial distributions of the IKN GCs 


Investigating whether the mass of the IKN GCs varies as a func¬ 
tion of the galactocentric distance may hold important informa¬ 
tion about cluster formation and migration due to dynamical fric¬ 
tion. Here we attempt to derive a new center of IKN from the 
stellar positions of its RGB stars. 

Deriving the radial distribution of the GCs from the IKN 
centre is a rather complicated task because of the partial cov¬ 
erage of the HST/ACS field and the strong incompleteness due 
to the bright foreground star. To derive the photometric center 
of IKN, we used the positions of the field stars detected in the 
HST/ACS image whose photometric and spatial data was pro¬ 
vided to us (by jLianou et al.|2010| private communication). In 
Figure |7] we show the two-dimensional number density distri¬ 
bution of the IKN field stars. To minimize the impact of the 
foreground star in deriving the IKN photometric centre, our ap¬ 
proach is first to obtain the histogram distributions along the x- 
and y-axes of the ACS field and then to fit a Gaussian functions 
to each. Their peaks positions are adopted as the IKN photo¬ 
metric centre. The IKN photometric centre derived using this 
method is shown with a white star in Figure|7] with coordinates 
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Fig. 8 Globular cluster luminosity versus projected distance form 
galaxy photometric center for GCs in dwarf galaxies (small grey 
circles) and IKN known GCs (large blue squares) and IKN GC 
candidates (large open squares). The running median (red dia¬ 
monds) is computed from the entire dwarf galaxies GC sample 
with a bin size of 0.345 Kpc, where the y-errorbar is the standard 
deviation of the number of GCs in each bin. 


of RA = 10 : 08 ; 05.8 and DEC = -h 68 : 25 : 20 and with a 
total uncertainty of 8.5"(150pc) derived from the fit. We note 
that the probably the choice of a Gaussian function may not 
be optimal, but given all other involved uncertainties, it is suf¬ 
ficient for our analysis. Additionally, to obtain a more realistic 
constraint to the IKN center we fit ellipses to the arcs defined 
by the constant surface stellar density contours. The 2D surface 
density of the RGB stars was estimated in a running window 
with a width of 300pixels (15")- The fit was performed for the 
arcs that are 300 pixels inwards from the ACS detector edges 
and with constant surface density between 3 and 6 stars /arc- 
sec^. The results is that the center is at RA= 10 : 08 : 24 + 45" ; 
DEC= -h 68 : 24 ; 47 + 73" , PA= 154.4 +11.1 degrees (position 
angle east of north), ellipticity= 0.216 ± 0.12 (e = 1 - bja). 


We can use center of IKN to see how the GCs are distributed 
as a function of the projected galactocentric distance. The ra¬ 
dial distribution of the IKN GCs is shown in Eigure[^with solid 
(dark blue) squares, where we also show the projected radial dis¬ 
tributions of GCs from a larger sample of dwarf galaxies from 
[Georgiev et al.| ( |2()08[|2009^ . In general in Eigure[^is seen that 
the cluster luminosity (mass) decreases with increasing galacto¬ 
centric distance. The high mass clusters are closer to the center. 
To quantify this trend for the combined Rpjoj distributions from 
all dwarf galaxies, we calculate the running median as a func¬ 
tion of Rproj with a bin size of 0.345 kpc (diamonds connected 
with a solid line). Here we confirm an earlier finding that more 
luminous clusters in dwarf galaxies are typically more centrally 
located ( |Georgiev et al.|2009a| ). Clearly, in Eigurej^is seen that 
the most massive IKN GC is closest to the derived here pho¬ 
tometric center. A likely explanation is that dynamical friction 
may have “sorted” the IKN GCs. However, spectroscopic obser¬ 
vations are required to know the local velocity field and whether 
dynamical friction can be an efficient process for IKN central 
stellar density. 
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3.5. Alignment of IKN GCs and orientation in the M81 group 

We also report that the IKN GCs projected positions appear to 
be linearly aligned. Although this could be also a chance align¬ 
ment due to low number statistics, the exact spatial and kine¬ 
matic structure of the IKN GCs needs to be spectroscopically 
confirmed, because it can have an important implication for in¬ 
terpreting the formation origin of IKN, its SEH and interaction 
within the M81 group of galaxies. Eurthermore, a similar lin¬ 
ear alignment has been previously reported for the old GCs in 
NGC 6822 ( |Hwang et al.|20TT]|Huxor et al.|2013| l. 

The projected linear alignment between the IKN GCs can 
be already seen in Eigure[T] In Eigurej^a) we show the straight 
line least squares fit through the GCs’ positions. Eor reference, 
we show with an ellipse the size of IKN’s Holmberg diame¬ 
ter, 2.69 kpc, (measured at pg - 25 mag) and its ellipticity, 
e - I - bja - 0.18, retrieved from the HyperLEDA and NED 
databases. For the center of the ellipse we used the center we ob¬ 
tained in § 3.4| We note however, that the relative uncertainty in 
the IKN ellipticity and particularly in the orientation of its major 
axis reaches nearly 100%. This thus prevents us from a more re¬ 
liable discussion on the alignment between the IKN major axis 
and the fitted line through the positions of its GCs. 

In the bottom panel of Figure.j^we compare the orientation 
of the GCs alignment plane with respect the large scale distribu¬ 
tion of galaxies within the M81 group. If the GCs alignment is 
due to their formation in a disk or triggered by interactions with 
galaxies in the M 81 group, angular momentum conservation ar¬ 
guments would suggest a correlation with the positions of the 
nearest galaxies and the larger structure in the group. Figure]^ 
shows that the GC’s projected alignment and the IKN major axis 
seem to be orthogonal to the larger scale galaxy distribution, i.e. 
the plane connecting IC 2547-NGC 3077-M 81. However, spec¬ 
troscopic observations are required in order to properly describe 
the kinematics of the IKN GCs, whether they are indeed on a 
(rotating?) disk and what is the orientation of the vector of the 
angular momentum. 


3.6. Implications for the formation origin of IKN 


Hints of the unusual and probably complex star formation his¬ 
tory of IKN were first highlighted by the presence of a large 
number of old globular clusters compared to its total luminos¬ 
ity, i.e. the highest known GC specific frequency of = 126 
( |Georgiev et al.||2009b| |2010[ ). Subsequently, resolved stellar 
population HST/ACS CMD analysis ( [Lianou et al. 2010 1 , re¬ 
vealed that IKN has a very broad metallicity distribution func¬ 
tion. It peaks at a much higher metallicity for the IKN lumi¬ 
nosity, unlike other early-type dwarf galaxies in the M 81 group 
with similar luminosity ( |Lianou et al.|2010) l. This led those au¬ 
thors to consider IKN as a galaxy which may have formed as a 
tidal dwarf galaxy during galaxy-galaxy interactions in the M 81 
group. Recently, Larsen et al. ( 2014a| l performed a direct spec¬ 
troscopic comparison between the brightest IKN GC-5 and those 
of another faint dwarf galaxy - WLM, who showed that in both 
galaxies the GCs are metal poor, although the IKN-05 GC spec¬ 
tra is of a very low S /N. 

Our analysis here further strengthens the idea of a complex 
SFH for IKN, based on its GCs’ ages, metallicity and spatial dis¬ 
tributions. Although the photometric uncertainties render large 
errors in these quantities, the photometric ages and metallici- 
ties of the GCs places them among the metal poor field stellar 
component (cf Fig.[^. This suggests that the GCs and the metal 
poor field stellar component have formed during the same burst 
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Fig. 9 Top: Spatial distribution of the IKN GCs (large solid cir¬ 
cles). The htted linear regression through their position is shown 
with solid line, while the IKN major axis and orientation are 
shown with dashed and light gray line and ellipse, respectively. 
IKN RGB stars are shown with small (red) dots. Bottom: Po¬ 
sitions of galaxies associated to the M81 group according to 
Karachentsev et al. ( 2002| l. Circle sizes are proportional to the 
galaxy absolute B-band luminosity. The line passing through 
IKN roughly shows the orientation of its GCs from the top panel. 


of star formation. Our analysis of the field RGB star MDF sug¬ 
gests that there must have been at least two major star forma¬ 
tion episodes in the IKN SFH. If galaxy-galaxy interactions in 
the M 81 group have triggered bursts of intense star formation, 
likely galactic pericenter passages were responsible for the mul¬ 
tiple peaked SFH. If IKN has formed as a tidal dwarf galaxy, one 
would also expect that it has a higher for its luminosity metallic- 
ity, which is the case. The apparent alignment of the IKN GCs, if 


real, also would suggest that IKN could be a dSph transformed 
from a late-type disk dwarf. 

The relevant observational evidence for a complex SFH of 
IKN are; i) on average slightly more metal-rich GCs [Fe/H] = 
-1.4dex (cf. Fig.[^, than typical for a truly metal-poor GC pop¬ 
ulation of [Fe/H] « -1.7 dex, e.g. in the Fornax dSph all GCs 
are with [Fe/H] < -2.0 dex ( [Larsen et al.|2014bl [Hendricks et al.| 
[2014] l, ii) non-unimodal MDF of the field RGB stars combined 
with the apparent alignment of the IKN GCs place constraints 
to the possible scenarios that could have led to the complex IKN 
SFH. Namely, 1) the IKN GCs should have formed together with 
the main galaxy stellar population, so that the GCs’ and the RGB 
metallicities and MDFs are matched (Fig.[^; 2) this must have 
happened more than 5-8Gyr ago; 3) the intensity of star forma¬ 
tion must be such that it leads to the formation of an unusually 
large number of GCs or 4) IKN must have lost at least 50% of its 
initial mass in order to lead to the high GC specific frequency. 
Tidal formation origin and/or the tidal interaction triggered star 
formation scenarios can match these observations. However, if 
IKN is a primordial dwarf it must contain a significant amount 
of dark matter, which cannot be estimated for the currently avail¬ 
able data. Clearly, a spectroscopic follow up is needed to assess 
the dynamics of the system. 


4. Summary and conclusions 


We present a new K-band photometric analysis of the globular 
clusters of the IKN dwarf spheroidal using data taken with the 
NIR detector (LIRIS) installed on the 4.2m William Herschel 
Telescope. Combined with existing HST/ACS V and I data, we 
effectively resolve the age-metallicity degeneracy in the /-K^ vs 
V - I colour space. The new Kg -band photometry is tabulated 
in TablejH 

• Using SDSS DRIO archival images we performed a PSF 
photometric analysis to look for GC candidates in the field 
around IKN that is not covered by the HST/ACS observations. 
We found two new GC candidates based on departure form a 
stellar PSF (cf Fig.[^ § [2.2.3 1. One of the candidates (GC-6) is 
however very close to a bright star which renders its measured 
parameters highly inaccurate. High spatial and spectral resolu¬ 
tion and multi-wavelength observations are required to reveal the 
nature of these two objects. If one of these GC candidates proves 
to be another genuine GC, this will make IKN the galaxy which 
formed GCs with a very high efficiency per unit luminosity with 
one of the highest specific frequencies known to date, S n - 126. 

• To derive GCs’ photometric ages, masses and metallicities 
for all known GCs in the IKN galaxy in §[3.1[we interpolate be¬ 


tween Bruzual & Chariot ( [2003| l SSP model tracks (cf Fig.|^. 
For all GCs we derive mean photometric age which is consistent 
with the oldest ages in the SSP model grid (14.29/[ Gyr, Fig.|^ 
Table|^. The old GC age indicates an ancient burst of star for¬ 
mation. Using the empirically established observational relation 
between the most massive cluster and the galaxy star formation 
rate (see § [3.2| l, we find that the burst of star formation led to the 
formation of 3 x 10 ^Mq GC was about 10 M^lyear at about 14 
Gyr ago (cf Fig.|^ in §[3[2i. 

• We derive GC metallicities in the range 2.0 < [Fe/H] < 
-1.2 dex (Tab.|^ § j3.3| l. These are consistent with the metal-poor 
component of the metallicity distribution function of the field 
RGB stellar population of IKN reported by Lianou et al. ( 2010[ l. 
Our reanalysis of the field RGB stars MDF, data kindly provided 
to us by Lianou et al.| confirms their finding of a broad MDF. 


However, we find at least two major bursts of star formation 
with a small statistical significance, the most recent of which 
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produced the most metal rich stellar population with a peak of 
[Fe/H] -l.Odex (Fig.|^. Although the photometric metallic- 
ities of the IKN GCs are large, they are more consistent with 
the metal-poor component. Our mixture model analysis suggests 
that about 72% of the IKN stellar population could have formed 
at the same star-formation epoch as the bulk of the IKN GCs. 
We note however, that a spectroscopic confirmation is necessary 
in order to firmly establish the IKN metallicities. For GC-5 such 
a confirmation does exist and it falls within the photometric un¬ 
certainty estimate. 

• We find that the more luminous (massive) IKN GCs are 
preferentially distributed toward the center of IKN (§ 3.4 1 , con¬ 
sistent with earlier observations (cf Fig.[^. Such mass segrega¬ 
tion has been suggested to arise due to dynamical friction argu¬ 
ments, which is efficient in low velocity systems, but strongly 
depends on the local stellar density. Further spectroscopic obser¬ 
vations are required to map the velocity field of IKN and estab¬ 
lish whether dynamical friction wold be efficient for such a low 
stellar surface density galaxy. 

• We find that the projected positions of the IKN GCs appear 
to be aligned. (Fig.[^ §|3.5[). This is similar to the observed align¬ 


ment of the GCs in the Local Group dSph NGC 6822 (Hwang 
|et al.|[20TT| |Huxor et al.||2013'] l. It remains to be shown spec¬ 
troscopically whether the IKN GCs have disk kinematics or the 
apparent distribution is a chance alignment. Knowing this will 
help for the interpretation of the IKN origin, as a merger of dwarf 
galaxies will destroy existing disk structures ( Bekki|2008) l. The 
orientation of the plane of alignment of the IKN GCs appears to 
be very close to the orientation of the IKN major axis, however, 
the latter is very poorly constrained. Looking at orientation of the 
GC alignment compared to the large scale distribution of galax¬ 
ies reveals tentatively that it is orthogonal to the large plane con¬ 
necting several of the M81 group galaxies (see § |3.5| and Fig.[^ 
bottom panel). 

In summary, the large cluster formation efficiency per galaxy 
mass, high GC metallicities and positional alignment, point to¬ 
wards a complex SFH for IKN. Whether the tidal origin scenario 
for the galaxy indeed holds remains to be cleared through further 
spectroscopic investigations. 
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Table 4. Photometric properties of the IKN globular clusters (GCs) and the new candidates (GCCs). 
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